A Solution for Faults at Two Locations in Three-Phase 

Power Systems 

By E. F. VAAGE 

This paper is an outgrowth of studies of double faults to ground 
in three-phase power systems made by the author in connection 
with work of the Joint Subcommittee on Development and Re- 
search, Edison Electric Institute and Bell Telephone System. The 
paper provides a systematic solution, based on the method of 
symmetrical components, by means of which currents and voltages 
can be determined at times of fault involving any combination of 
phases at one or two locations on three-phase power systems. 

1. Introduction 

A KNOWLEDGE of the magnitude and phase relation of power 
system voltages and currents for various types of faults in three- 
phase systems is of importance in the study of various problems, among 
which are relaying studies, the efhcacy of current limiting devices and 
their reaction on the power network, and estimates of induction in 
paralleling communication circuits. 

The method of symmetrical components as developed by Fortescue ^ 
and others is now extensively used in the solution for currents and 
voltages in three-phase power systems under fault (short-circuit) 
conditions. Formulas for special cases of faults, such as single and 
double line-to-ground faults, can be found in various text books on this 
subject. The solution for simultaneous faults at two locations has 
been treated by Miss Clark, ^ in a form particularly adaptable to the 
use of a calculating board. 

The present development provides a complete and systematic solu- 
tion for currents and voltages at times of fault on any number of 
phases at one or two locations in a three-phase system, in which gener- 
ators may be assumed in phase and of the same internal voltage, and 
where load currents can be neglected. These are the usual assumptions 
mads in computing fault currents, except for certain special problems, 
such as that of power system stability. The methods employed herein 
could be extended to cases where generators of different phase angles 
and voltages of more than two points of fault are involved. Formally 
such cases can be treated in a manner similar to that given in the 
paper. The number of impedances to which an n-terminal network 

• Reference numbers refer to references appearing at the end of the article. 
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can be reduced is given by the expression \n{n — 1). For « = 3, 
the case treated in this paper, three impedances are required which 
necessitate six equations for the general solution of the six fault 
currents. For » = 4, which would be the case for three points of 
fault (or two points of fault and two generating voltages), six im- 
pedances would appear in the reduced network and this would necessi- 
tate twelve equations for the general solution of the fault currents. 
For larger values of n, the necessary number of equations increases 
rapidly, thus making the solution impractical. Such problems 
usually, as a practical matter, are more readily solved by the use of a-c. 
calculating boards. 

While no departure from the general methods of symmetrical com- 
ponents is made in the present development, a systematic method of 
handling the equations is presented and means of determining the 
coefficients given so that numerical calculations can be directly carried 
out when the constants of the network are known. 

2. General Solution 
The equations developed in this paper are based on the sequence 
impedances looking into a three-phase network from two points of fault. 
Consider the network shown in Fig. 1. This system can be reduced 



CONTArNS ENTIRE NETWORK 





FAULT A FAULT B 

Fig. 1 — General network diagram. 

to an equivalent star for each of the positive, negative and zero se- 
quence networks, with legs to the points of generation and to the faults 
at A and B. Figure 2 shows the reduced positive sequence network. 
Similar diagrams can be made for the negative and zero sequence 
systems except for the fact that in these cases there are no generated 



292 



BELL SYSTEM TECHNICAL JOURNAL 



voltages, and the impedances and currents are the negative and zero 
sequence quantities. 

The reduction of a network to an equivalent star is usually a tedious 
and sometimes a difficult process especially in large interconnected 
systems. Methods of accomplishing the reduction, such as delta-star 
transformations, simultaneous equations or direct measurements on 
calculating boards can be found in the literature.^' ^ 

Having reduced the three sequence networks to equivalent stars, 
the equations are developed as shown in the Appendix. 

The following set of equations (1) is the general solution for fault 
currents during simultaneous three-phase faults to ground at two 
different locations in a power network. The set applies directly to 
the calculation of ground fault currents on a system having finite 
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fau'ltA faults 

' Fig. 2 — Reduced positive sequence diagram. 

neutral impedances or an isolated system in which the zero sequence 
capacitance is taken into account. For other types of faults, such as 
faults to ground in isolated systems in which zero sequence capacitance 
has been neglected, or for phase to phase faults, set (1) is not directly 
applicable since some of the constants become infinitely large. How- 
ever, by certain transformations of set (1), more convenient sets (2) 
and (3) are obtained, directly applicable for solution of these latter 

cases. 
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The six equations are written In matrix form with the currents and 
voltages outside the system matrix. For example the first row in (1) 
is interpreted as : 

AiJa^, + -4 12/^6 + AiJac + Aulsa + Aislsb + AuIbc = SE 

The values of the A 's in (1) are given in Table I. It should be noted 
that of the 36 constants only 13 are distinct. Six of these are in the 
nature of self-impedances, two are transfer impedances between phases 
at A and two between phases at B. The remaining three are transfer 
impedances between the two faults at A and at B. 

Considerable reductions in the constants are obtained when the 
positive and negative sequence impedances are assumed equal. These 
values are given in Table 11. 

Faults to ground on less than three phases at one or both locations 
are accounted for by assuming the corresponding fault resistances in- 
finitely large. The currents to ground in the sound phases are zero. 
Striking out the columns containing these currents and the correspond- 
ing rows, indicated by the index at right in equation (1), a reduced 
set of equations is obtained from which the desired currents can be 
found. A few examples are given in subsequent sections. 

In power networks with isolated neutral the zero sequence impedance 
Zco reduces essentially to the capacitance of the system. In this 
case equations (1) are still appropriate and will give a rigorous solution 
for the six currents. However, in many cases it is sufficiently accurate 
to neglect the capacitance of the system. This results in infinitely 
large values of all of the A's in Table I (Table II), since each depends 
on Zco which is infinitely large. For this condition it is desirable to 
transform the set of equations in (1) to a more convenient set with 
finite constants. 

The transformation required is obtained by observing that, with 
Zco infinitely large, the sum of the zero sequence currents Iao + /bo 
must be equal to zero. Making use of this relation the difference of 
the zero sequence voltages at A and B (equation (50) of appendix) 
reduces to : 

Vaq — Vbo = {Zao + Zbo)Ibo 

The last equation shows that subtraction of equations associated 
with phases at A from those at B removes the infinitely large element 
Zco. This can be done in nine ways (ignoring reversals of sign), but 
three of these result in the single equation : 

Iah + IaI + I Ac + Isa + iBb + I Be = 

This equation with any five of the remaining six constitutes an inde- 
pendent set; for convenience in dealing with special cases the redundant 
set of seven equations is shown in the following array : 
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Isolated System — Capacitance Neglected 
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(2) 



I A. + lAb + Iac + Tea + /b5 + /b. = (2a) 

The index to the right indicates which of the equations in (1) have 
been used. The values of the B's are given in Table I and Table II. 
In case of faults to ground on less than three phases, as in equations 
(1), columns and rows associated with sound phase currents are to be 
deleted; with respect to the rows, however, the index is double and 
all rows having the index of the sound phase or phases are deleted. 
If, for example, the sound phase is Aa, rows 1 and 2, each of which 
contains Aa in its index, as well as column 1, are deleted. This leaves 
only four equations, which together with (2a) give the necessary five 
equations for the five currents. For this reason all six equations are 
given in (2), since any phase might be involved in special cases. 

Phase-to-phase faults are obtained from the general case (1) by 
allowing the resistances Raf and Rbf to become infinite. In this case 
phase-to-phase quantities at the same location remain finite and the 
appropriate set of equations is obtained by subtracting equations 
having the corresponding phase indexes; thus Aa — Ah, Aa — Ac and 
Ah — Ac indicate subtractions at A. There are six possible ways of 
doing this, ignoring reversals of sign. The resulting set is given in (3). 
The four equations obtained by taking any two of the first three and 
any two of the last three equations in this set together with the two 
equations (3a) relating to the sum of the currents at each fault location, 
which from physical considerations equal zero, constitute an inde- 
pendent set. For convenience in dealing with special cases all eight 
equations are given below: 

Phase-lo-Phase Faults 
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The index to the right indicates which equations of (1) have been 
used. The values of the Cs are given in Table I and Table 11. 
The total ground fault currents at the two fault locations are : 

IaF = lAa + I Ah + Ia^ (4) 

/bF = Isa + /b6 + Ibc (S) 

and the total residual current in the two faults is; 

Ib = Iaf-\- Ibf (6) 

In an isolated system in which capacitance has been neglected this 
current is zero and equation (6) will be identical with (2a). 

The distribution of these currents in the network can be found as 
follows. From equations (51) and (52) of the appendix the calculated 
fault currents are transformed into sequence currents. By working 
back into the original sequence networks the sequence currents in each 
branch of the system'can be found and later combined from similar 
expressions as shown in (43) and (44) to obtain the actual branch 
currents. ^ 

A combination of the equations in (1) and (3) can be used for cases 
involving faults to ground at one location and faults between phases 
(not involving ground) at the other location. 

For faults to ground at A and between phases at B, the three first 
equations in (1) and the three last in (3) together with the last in (3a) 
constitute the most convenient set of equations for this type of fault. 

It should be noted that if all three phases are involved at B any two 
of the three last equations in (3) together with the last in (3a) can be 
used, while for less than three phases involved, the rules for striking 
out rows and columns automatically will result in the proper equations 
to be used. 

Vice versa the three first equations in (3) together with the first in 
{Za) and the three last equations in (1) will give the solutions for 
phase-to-phase faults at A and ground faults at B. 

This will be illustrated with an example in a later section. 

The voltages to ground at the two locations of faults can be obtained 
directly from equations (53) and (54) of the Appendix, after the 
currents have been evaluated. At any other point in the system the 
voltages to ground are found by adding the voltage drops of the lines 
in question to these voltages, treating each sequence network sepa- 
rately, then adding the sequence voltages together according to 
equations (46) or (47). 
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3. Special Cases 

The application of the three sets of equations (1), (2) and (3), will 
be illustrated with a few examples. For simple cases, such as a single 
or double line-to-ground fault at one location, the equations reduce to 
formulas frequently found in the literature on this subject. 

From set (1) equations for faults to ground at one or two locations 
can be obtained directly when the zero sequence impedance is finite. 
Set (2), obtained from (1), is the most convenient set for solutions of 
faults to ground in isolated systems in which capacitance has been 
neglected. The phase-to-phase fault currents are best obtained from 
set (3). 

3.1 Single Line-to-Ground Fault at A 

Consider a fault to ground on phase " 6 " at ^ . The solution can be 
obtained from (1) by letting: 

RAa = Rac = Rsa = Rsb = Rbo = ^ (7) 

This results in : 

Ia. = I Ac = Ib. = Isb = Ib. = (8) 

Striking out all columns in (1) containing the currents in (8) and 
the corresponding rows indexed by Aa, Ac, Ba, Bh and Be only one 
equation is left : 

^22/^6 = ^a-'E (9) 

The numerical value of A21 can be calculated directly from Table I, 
or on substituting the symbolic value of A^i in equation (9) the 
result will be : 

where 

Zi = Zai -\- Zci 

Z2 — Zaz + Zc2 ,. .^ 

Zq = Zao 1 Zco 
Rp = Rai + Raf 

Equation (10) is the well-known formula for a single line-to-ground 
fault at one location in a three-phase system. 

3.2 Double Line-to-Ground Fault at A 

Consider a double line-to-ground fault on phases "a" and "b" 
at A . Then : 

Rao = Raa = RBb = Rb.= ^ (12) 
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and 

Iac = Ib. = Isb = Ib. = (13) 

Striking out the columns of (1) containing the currents in (13) and 

the corresponding rows (Ac, Ba, Bb and Be) the following two equations 

remain : 

AiJau + AuIai = 3£ 

(14) 
A.Jau -\- A22lAb = 3am 

from which on substituting the numerical values for the A's from 
Table I the two currents Iao and Jai can be found. The total fault 
current to ground at ^4 is: 

IaF = lAa +lAb (15) 

In the special case where Rao, R-Ab and Raf are zero the expression 
for Iaf can be reduced to the following expression after a direct sub- 
stitution for the ^'s in (14) is made : 

7 = - 3aZ2-£ ,jg. 

Z(iZ\ -\- ZqZ2 ~r ^1^2 
where : 

2i = Zai + Zci 

Z2 = Za2 + Zc2 (17) 

Zo = Zaq + ZcQ 

3.3 Simultaneous Double Line-to-Grotmd Fault at A and Double 
Line-to-Ground Fault at B 
Consider a fault-to-ground on phases "a" and "b" at A and phases 
"a" and "c" atB. Then: 

Rac = i?M = «> (18) 

Hence : 

Iac = lBt = (19) 

Striking out the two columns containing Iac and Isb and the two 
corresponding rows (Ac and Bb), the four following equations remain : 

AniAa + AizIai + AiJbo + AigIbc = 3£ 

AziIao. + A22lAb + A2JBa + Az^Bc = 30^£ 
AiilAa + ^42-^^/16 + A ulna "f- Aiz^Bc = 3£ 

AiiilAa + -462-^^(1 + ^fii-^Bu -f- -^es-ffic = 3a£ 
A symbolic solution in terms of the sequence impedances for these 
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currents becomes quite involved and it is advisable to substitute 
numerical values of the constants before solving for the four currents. 
The total fault currents at A and B, respectively, are (from (4) and 
(5) in connection with (19)) : 

Iaf = Ia. + lAb (21) 

Ibp = iBa + Tbc (22) 

In a similar manner faults to ground for any other combination of 
faulted phases can be found. 

3.4 Single Line-io-Ground Faults at A and B in an Isolated System 
Consider a fault-to-ground on phases "a" at A and "h" at B in 
an isolated system in which capacity can be neglected. Then : 

RAb = Rac = RBa = Rbc= <=^ (23) 

and 

Iaj. = Iac = Ib. = Ibc = (24) 

Striking out the columns of (2) containing the currents in (24) and 
the corresponding rows Aa — Be, Ab — Ba, Ab — Be, Ac — Ba. and 
Ac — Bb (all rows containing Ab, Ac, Ba and Be), leaves only one 
equation in (2), which together with {2a) gives: 



BiJa. + BiJbi> = 3(1 - a')E 
lAa -\-lBb = 

Solving for these currents the result is : 



(25) 



-On — £>u 
Inserting the values of the B's from Table I this reduces to : 

Zu = Zai -\- Zbi + 3Zci 
Z2i = Za2 -\- Zb\ + 3Zc2 
Zoi = Zao + Zb^ (28) 

Ra = RAa 4- Raf 
Rb ~ Rsb + RbF 

The subscript i (isolated) is used to distinguish these impedances for 
the isolated system from those used in (11), (17) and (38). 
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3.5 Phase-to-Phase Fault at A 
Consider a fault between phases "a" and "b" at A. Then let: 

Iac = iBa = Ish = Ibc = (29) 

Striking out the columns in (3) containing the currents in (29) and 
the corresponding rows (all rows containing Ac, Ba, Bb and Be) only 
one equation is left : 

CxilAa + CnlAt = 3(1 - a=)£ (30) 

It is further known from {J,a) that : 

lAh = - lAa (31) 

Substituting (31) and the constants Cn and Cu from Table I in 
(30) the result is: 

iAa - lA^ - 2j + Z, + Ra. + Ra, ^^^^ 

Z\ = Za\ + Zc\ .,„ 

Z2 — Za2 -\- Zci 

which is a well-known expression for a phase-to-phase fault. 

3.6 Three-Phase Fault at A 
For this case : 

Ib. = Isb = Ibc = (34) 

Striking out the columns in (3) containing the currents in (34) and 
the corresponding rows, three equations remain, any two of which 
together with the equation from (3a) relating to the currents at A give : 

CuIac + CnlAb + CuIac = 3(1 - a^)^ 

Czi/xa + C22/AJ. + CislA. = 3(1 - a)E (35) 

lAa + lAb + Iac=0 

from which the currents can be found. 

In the special case where the fault resistances are all zero, the three 
currents are equal in magnitude and related as follows : 

lAa = alAh = a^Ac (36) 

The rank of the system determinant in (35) is therefore 1. Using 
any of the first two equations in (35) in connection with (36) and the 
constants in Table I, the result is: 

lAa = alAb = aUAc = -™- (37) 

Zr = ZAii-Zci (38) 



302 BELL SYSTEM TECHNICAL JOURNAL 

3.7 Phase-to-Phase Fault at A and PJiase-to-Ground Fault at B 
Consider a fault between phases "a" and "b" at A and a fault to 
ground on phase "c" at B. Then : 

Ra. = Rsa = Rsb^ '^ (39) 

and , , 

Ia. = iBa = lnb = (40) 

As explained in a preceding section the three first equations in (3) 
together with the first in (3a) and the three last equations in (1) may 
be used for this case. 

Striking out the columns Iac, Ib^ and Isb and the corresponding rows 
Aa — Ac, and Ah — Ac in the three first equations in (3) leaves only 
the first equation. Similarly by striking out the columns Isa, iBb, and 
the corresponding rows Ba and Bb in the three last equations in (1) 
leaves only the last equation. Hence: 

CuIag + CulAb + CibIbc = 3(1 — a')E 
Aei^Aa + A62lAb + Aaelsc = 3a£ 

and finally from (3o) : 

lAa + lAb=0 (42) 

from which the three currents can be found. The A's and C's are 
given in Table I and Table II. 

4. Conclusion 

While the probability of all phases being faulted at both locations 
simultaneously is very remote, the three sets of equations (1), (2) and 
(3) have been given in such a form that they conveniently will provide 
a solution for any combination of phases faulted from a single line-to- 
ground fault at one location to the most involved fault condition. 

In Section (3) of this paper, in which special cases have been treated, 
only simple types of fault conditions have been shown in order to 
illustrate the method to be used and to prove that the general equations 
reduce to well-known formulas. 

The constants given in Table I consist of the nine quantities S^, 
Sb, So, Ta, Tb, Tc, Ua, Ub and Uc arranged as shown for each set of 
equations. Table II gives somewhat simpler values for the constants 
in cases where the positive and negative sequence impedances are 
assumed equal. 

The voltages to ground at the two fault locations are given by (46) 
and (47) in the Appendix. 

It is hoped that this development will provide a more unified presen- 
tation of fault current calculations in power networks. 
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APPENDIX 

The standard notation for phase and sequence quantities is usually 
indicated by a subscript. Thus la, lb, etc., means the current at the 
point of fault of phase "a" and "6," respectively, /i, a^Ii and ali 
are the positive sequence currents in phase "a," "b" and "c," re- 
spectively. In this treatment, however, complication arises from the 
fact that two points of faults are involved and it will be necessary to 
distinguish between the quantities at these two locations. This is 
most conveniently done by a double subscript, the first referring to the 
point of fault and the second to the phase or sequence in question. 
Thus Iag, Isa, etc., are the currents in the fault at A and B of phase 
"a" and Iai, Ib\ the positive sequence current at the two points of 
fault, respectively. Making use of this notation the fundamental 
equations for the sequence currents at fault A are : 

■^xo + Iai + Iai = Iao 

Iao + a'^lAi + alA-i = Iai (43) 

7^0 + uIai + c^/a2 = Iac 
And at fault B : 

Ibo ~I~ Ibi + Ib2 = Isa 

Ibo + a^Ini + a/fla = Isb (44) 

Ibo + alsi + a^/ea = Ibc 

where the coefficient "a" is the sequence operator, having the value: 

1 , .V3 

1 .V3 (45> 



'^=-2+-?" 2 



"'=-2--^T 

The voltages to ground at the two fault locations are given by : 

VAa = Vao + Vai + Va2 

— iRAa + -^Ap)/aq + RAFlAb + RafI Ac 

VAb= VA. + a'^VAi + aVAi. 

= RafIa. + {RaU + RAF)lAb + RafIac ^ ^ 

Vac = VAo + aVAi-\-a?VA2 

= RafIa^ + RAFlAb + (Rac + Raf)Iac 
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VBa = Ffio + Vbx + Vb2 

= iS-Ba + RbFJIbb H~ RsplBb + RbfIbc 

Fflfc = Fbo + a^Kfii + aVBa /47J 

Vbc == FBo + aFBi + a''l^B2 

= RbfIbh + RBF^Bb + (-Kflc + RBF)^Be 

Consider the positive sequence diagram in Fig. 2. Evidently: 

Vai = E — {Zai + Zci)!ai — ZciIbi ,^^\ 

Vbi = E — ZciIai — {Zb\ + Zc\)Ib\^ 

where Va\ and F^bi are the positive sequence voltages to ground at the 
two fault locations. Similar expressions can be obtained for Vau, 
Vb2, V^o and Vbq, except for the fact the E is zero in these cases and 
the impedances and currents are the negative and zero sequence 
quantities. They are : 

Vaz = — {Zaz + Zc2)Ia2 — ZcilBt ,^g, 

1^52 = ~ Zc2^A2 ~ {Zb2 + Zc2)!b2 

Vao = — {Zao + Zco)!ao — ZcoIbo ,-Q^ 

Vbq = ~ Zco^Ao — {Zbo + ZcojIbq 

Solving (43) and (44) for the sequence currents the result is : 

Iao = \{lAa + ^Ab + Iac) 

Iai = 1(Ia. + alAb + a'lA.) (51) 

Ia2 = ^(^Aa + a'^lAb + alAc) 
IbO = si^Ba + ^Bb + Ibc) 

Ibx = Wb. + aiBb + aHBc) (52) 

Ib2 = Wb. + aHBb + c/b.) 

Substituting the expressions for the sequence currents in (48), (49) 
and (SO) the result is: 

Vai = E- \{Zai + Zci){Ia. + alAb + a'/x.) 

- \Zci{lBa. + aiBb + oVbc) 
Fa2 = - 1{Za2 + Zc2)(I^a + a^-^Aft + a/x.,) /e^-, 

~ aZco{-^Ba + /flb + -fjJc) 
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- K^Bi + Zci){Ib. + alsb + aHsc) 

Vni = — \Zci{lA<i + a^^Ab + alAc) fe^A\ 

-l{ZBi + Zc2){lBa^-anBb + alB.) ^ ' 

Vbq = — \Zc\i{lAa + I Ah + Iac) 

— K^BO + Zco){^Ba + i^B6 + ^Bc) 

Substituting (53) and (54) in (46) and (47) the six original equations 
in (1) are obtained. 

References 

1 C. L. Fortescue; "Method of Symmetrical Coordinates Applied to the Solution of 
Polyphase Networks," Trans. A. I. E. E., 1918, Vol. 37, pages 1027-1114. 

2. E. Clarke: "Simultaneous Faults on Three-Phase Systems," Trans, A. I. E. E., 

1931, Vol. 50, pages 919-941. 

3. Wagner and Evans: "Symmetrical Components" (McGraw-Hill, 1933). 

4. Monseth and Robinson: "Relay Systems— -Theory and Application" (McGraw- 

Hill, 1935). 
5 Lyon: "Applications of the Method of Symmetrical Components" (McGraw- 
Hill, 1937). 

6. G. Oberdorfer: "Der Erdschluss und seine Bekampfung" (Springer, 1930). 

7. M. J. Fallou: "Propagation des Courants de haute Frequence Polyphases le long 

des Lignes A6riennes de Transport d'Energie Affectees de Courts — Circuits 
ou de Defauts d'isolement," BtUletin de la Societe Frangaise des Eledridens, 
August, 1932, pages 787-864. 



